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level (CHPT)g with linear realization of chiral U(3) x U(3) symmetry. It is shown 
that for the explanation of the experimental probability of the — ► 7r + 7r + 7r~- 
nonresonant decay one does not need to use unnaturally heavy light current quarks 
as has been suggested in Ref.[l]. 
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1 Introduction 



The problem of the — > (7r + 7r + 7r~) NR -nonresonant (NR) decay has been recently dis- 
cussed in Ref.[l]. According to the conclusion suggested in Ref. [1] the theoretical expla- 
nation of the experimental probability of this decay can be reached via the adoption of 
unnaturally heavy light current quark masses m — \ {rn§ u + m od ) = 22 4- 38 MeV [1] 
instead of the widely accepted value m — \ (m 0u + m Qd ) = 5.5 MeV [2]. 

The application of the effective quark model with chiral symmetry involving Heavy 
quark effective theory (HQET) [3-5] supplemented by Chiral perturbation theory at the 
quark level (CHPT) g [6-8] with linear realization of chiral (7(3) x 17 (3) symmetry to the 
calculation of chiral corrections to the mass spectra and leptonic constants of charmed 
mesons [9], and to the form factors of semileptonic Z)-meson decays [10] has shown that 
the widely accepted values of the current quark masses m 0u = 4 MeV, m od = 7 MeV and 
rrios = 135 MeV [2,6-8] describe well the experimental data. 

This paper is to apply the effective quark model with chiral (7(3) x (7(3) symmetry 
incorporating HQET and (CHPT) g to the computation of the partial width and the 
probability of the Df — > (7r + 7r + 7r~)NR decay. Since experimentally the probability of this 
decay has been defined compared the probability of the 77+ — > 4>n + decay [11], i.e. 

we also compute the partial width of the 77+ — > 0tt + decay. 

The effective low-energy Lagrangian responsible for non-leptonic decays of the D+- 
meson reads [12] 

C cS (x) = V* V ud & [5/(1 - 7 5 ) c(x)} [« 7 „(1 - 7 5 ) d(x)\, (1.2) 

where Gf = 1-166 x 1CT 5 GeV -2 is the Fermi weak constant, V* s and V ud are the elements 
of the CKM-mixing matrix and C\ ~ 1.3 is the Wilson coefficient caused by the strong 
quark-gluon interactions at scales p > fi, where p, is a normalization scale. In (CHPT) 9 
we should identify p with the scale of spontaneous breaking of chiral symmetry (SB%S) 
A x = 0.94 GeV [7], i.e. p = A x = 0.94 GeV (see also [6,8]). The contribution of strong 
interactions at scales p < p = A x is described by (CHPT) g in terms of constituent quark 
loop diagrams, where the momenta of virtual quarks are restricted from above by the 
SBxS scale A x [6-10,13-16]. 

The effective quark model incorporating HQET and (CHPT) g resembles that suggested 
by Bardeen and Hill [17] that is also based on the Nambu-Jona-Lasinio model. There is 
only distinction that in the Bardeen-Hill model heavy-light mesons are considered like 
partners of light mesons, whereas in our effective quark model, heavy mesons are external 
states with respect to the light ones. This distinction influences only the redefinition of 
phenomenological parameters that are introduced in the model. Nevertheless, all results 
obtained within our effective quark model should be valid too in the Bardeen-Hill model. 
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2 The D+ 07T+ decay 

The amplitude of the non-leptonic decay — ► 07r + can be defined as follows 

M(£>+ - <M+) =< 7T+(g)0(Q)|£ eff (O)|L>+(p) > • (2.1) 

In (CHPT)g at the tree-meson approximation the computation of the matrix element of 
the non-leptonic decay Df — > <f)7t + Eq. (|2.1| ) agrees with the vacuum insertion approxi- 
mation leading to a factorized amplitude [7]. This gives 

M(D+ - 07T+) =< n + (q)<t>(Q)\/: eS (0)\Dt(p) >= 

= K Vud Cl < 7f+(9) 1 [U 7/1(1 ~ 10 > 

x <0(g)|[s 7M (l-7 5 )c(x)]| J D+(p) >= 

= — i Gp V; s V ud d g" < 0(Q)| [s^(l - 7 5 ) c(x)]|£>, + (p) > . (2.2) 

The computation of the matrix element < 7r + (g)|[u7 M (l — 7 5 ) d(x)]|0 > has been carried 
out in Ref. [7] 

< ir + {q)\[Wf(l - 7 5 ) d(a;)]|0 >= i V2F w q», (2.3) 

where F n ~ 92 MeV is a leptonic constant of the 7r + -meson [11]. 

The matrix element < 0(Q)|[s7 M (l — 7 5 ) c(a;)] |_D+ (p) > of the — > transition can 
be expressed in terms of the form factors of the semileptonic — > (piue decay [10,16] 

<<j>(Q)\[si,(l-l 5 )c(x)}\Dt(p) >= 

= iai(q 2 ) e;(Q) - % a 2 (q 2 ) (e*(Q) ■ p) (p + Q), - 1 a 3 (g 2 ) (e*(Q) • p) (p - Q)„ 

-2 6(9 2 ) (^° 123 = 1), (2-4) 

where e*(Q) is the polarization 4- vector of the 0-meson. As a result the amplitude of 
the .D+ — > 07r + decay reads 

M(L>+ -> 07T + ) = 

= G F K* s Kd Cx F w K(M^) - (M* + - Mj) o 2 M) - M 2 a 3 (M 2 )} q ■ e*(Q), (2.5) 

where we have set p 2 = M 2 D+1 Q 2 = M 2 and q 2 = M 2 . Setting M w = we reduce the 
r.h.s. of Eq. (|2.5| ) to the form 

M(D+ -»■ </>7T + ) = 

= G F Ci -Ftt [oi(0) - (M 2 Df - M|) a 2 (0)] g • e*(Q). (2.6) 

In HQET supplemented by (CHPT) 5 the form factors ai(0) (z = 1, 2) have been computed 
in Ref. [16] and read 

ai(0) = 1.47 GeV , a 2 (0) = 0.21 GeV" 1 . (2.7) 
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The partial width of the D~ s 



i )tt + decay is given by 



r(D+ - 07T+) = |c 1 | 2 |G F v;;K rf | 2 ^ 2 [a 1 (o) - (M 2 + - M 2 ) a 2 (0) 



2 5 



8vrM 2 



\C\\ 2 \G fV* s V ui ]\ 2 F% 

|Ci| 2 x 0.36 x 10 11 s -1 , 



ai (0)-(M 2 -M 2 )a 2 (0) 



M, 



2 \ 3 



64tt 



(2i 



where \q\ = (M^ + — M|)/2M<^ is a relative momentum of a 0-meson and a massless 
7r + -meson. The numerical value has been obtained at \V CS \ = 1.01, \V u d\ = 0.975, M D + = 
1.97 GeV, = 1.02 GeV [11]. We have dropped the dependence on the 7r-meson mass 
and set M w = 0. 

When using the partial width Eq. ( |2.8|) we can estimate r D + the time life of the D+- 
meson: 



B(D+ -> <t>^ 



(3.6 ±0.9) x 10 



-2 



D2 



|Ci| 2 x 0.36 x 10 11 (Cil 2 x 0.36 x 10 11 
(0.59 ±0.15) x 10~ 12 s, 



(2.9) 



where we have applied the experimental value B(D+ — ► </>7r + ) = (3.6 ± 0.9) % [11]. Our 
estimate r D + = (0.59 ± 0.15) x 10~ 12 s agrees well with the experimental value r D + = 
(0.467 ±0.017) x 10~ 12 s. 



3 The D+ (7r + 7r + 7r) NR decay 

The amplitude of the D+ — > (7r + 7r + 7r~) NR decay is defined 



< (7f + (q + )7f + (p + )7i-(p_)) NR \C cS (0)\Dj(p) >-- 
KV*dPi < (7r + (g + )7r + (p + )7r-(p_)) NR |w(0)7 M (l -7 5 )^(0)|0 > 



M(D+ -> (7T + 7T + 7r-) N R 

G p 
= ~V2 

x <0|s(0) 7 ' t (l-7 5 )c(0)|D+(p)> 



(3.1) 



The matrix element < 0|s(0)7 M (l — 7 5 ) c(0)\D£(p) > is determined by the leptonic 
constant F D +: 



< 0|s(0) 7 M (l -7 5 )c(0)| J D+(p) >= -iV2F Dt p», 

where F D + = 1.17 F D + = 1.41 F n [9,13]. Thus the amplitude of the D+ -> 
decay reads 

M(Di -> (tt+tt+Onr) = i G F V* V ud d F D + 

XJV < (7r + (g + )7r + (p + )7r-(p_)) NR |M(0)7 /t (l -7 5 )^(0)|0 > 



(3.2) 

(7r + 7r+7r~) NR 

(3.3) 



In (CHPT) 9 the matrix element < (7r+(g + )7r+(j9 + )7r-(p_)) NR |M(0)7 At (l - 7 5 )d(0)|0 > is 
defined by the box-constituent quark diagrams [8] and can be described by the momentum 
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integral 



< (7r + (g + )7r + (p + )7r-(p_)) NR |M(0)7^(l-7 5 )^(0)|0>=2-^ x 

16 TX Z 

d A k . { „ K 1 K 1 B 1 5 1 



■ tr<^ 7 M 7 5 7 5 s 7 5 - 
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7T 2 ^ [ ' ' m — k + p ' m — k + p — p + ' m — k + q + m — k 

+(<?+ ~ p+) + ■ ■ ■ , (3-4) 

where the ellipses denote the contribution of the matrix element of the vector current 
7/(0)7^(^(0) not contributing to the — ► (7r + 7r + 7r~)NR decay amplitude, then m = 
0.33 GeV is the constituent light quark mass calculated in the chiral limit [6-8], and 
g w = ^/2m/F n , the coupling constant of quark-pion interaction, satisfies the constraint 
Ng 2 /8ir 2 = 1 [13]. The former becomes rather obvoius, if one would take into account 
that F n = 0(y/~N) at iV — > oo and make a change F n — > F n yjN/3, where in the r.h.s. 
of this relation F n = 92 MeV. Then, one finds substituting the numerical data that 
Ng 2 /87c 2 = 3m 2 /4n 2 F 2 ~ 1. We should empasize that we have ignored the contibution of 
the 7r + -meson pole contribution to the r.h.s. of Eq. (|3.4|) . As it is shown in Appendix the 
contribution of the 7r + -meson pole is of order 0{1/ M / t+) compared with the contribution 

D s 

of the momentum integral left in Eq.( |3.4| ). 

Since in (CHPT) 9 [8] virtual momenta of constiuent quark loops are restricted by the 
SB^S scale, i.e. k < A x , and HQET assumes that M D + ~ M c , so in a heavy-quark mass 
limit we following the Appelquist-Carazzone theorem [22] reduce the r.h.s. of Eq. (|3.4j ) to 
the form 

< (n + (q + )n+(p + )n-(p^) NR \u(0) y(l - 7 5 ) d(0)|0 >= 

m 2 f d A k [ „ - 1 , 1 1 . . 

tr 7V - f - 7 5 f + (q+ <- p + ), (3.5) 



y/2F n M 2 D+ J 7r 2 i { m -k + q+ m-k 



where we have used the relations g n = y/2 m/F^, Ngl/Sit 2 = 1 and p 2 = M 2 + . The 
calculation of the integral over k gives one 

/ ^trj-yV ^-^Umg^^N. (3.6) 

J { m -k + q + m-k) Ng£ 



In the r.h.s. of Eq.( |3.6| ) we encounter a logarithmically divergent integral which is well de- 
fined in (CHPT)g by the compositeness condition g 2 I 2 (m) = 1 [6-9] (see also Refs. [25,26]), 
where we have used the definition of I^im) [6-9] 

h[m) - 16^7 rtK-Ff (3 - 7) 

and by applying the relation Ng 2 /%ii 2 = 1 [13]. The compositeness condition g^l 2 {rti) = 1 
realizes correct kinetic terms of low-lying mesons, gg-collective excitations, in effective 
chiral Lagrangians derived in effective quark models based on the Nambu-Jona-Lasinio 
approach [6-9,17,23,24]. 

Formally, the calculation of the momentum integrals representing one-loop constituent 
quark diagrams should be carried out keeping only the divergent parts and dropping the 
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contributions of the parts finite in the infinite limit of the cut-off. Such a prescription 
realizes a naive description of confinement of quarks. Indeed, dropping the finite parts 
of quark diagrams we remove the imaginary parts of them and suppress by this the ap- 
pearance of quarks in the intermediate states of low-energy hadron interactions. This 
naive description of confinement has turned out to be rather useful for the derivation of 
effective Chiral Lagrangians [6-9,23,24]. As has been shown in Ref. [21] this prescrip- 
tion can be justified in the multicolour QCD (N — > oo) with a linearly rising interquark 
confinement potential. Thereby, within such a naive approach to the quark confinement 
mechanism one can bridge quantatively the quark and the hadron level of the descrip- 
tion of strong low-energy interactions of hadrons. The cut-off A x = 0.94 GeV and the 
constituent quark mass m = 0.33 GeV can be considered in such an approach as input 
parameters and fixed at one-loop approximation via the experimental values of the pirn 
coupling constant g p and the leptonic pion constant F n [6]. Thus, we should accentuate 
that in the effective quark models based on the Nambu-Jona-Lasinio approach quark 
diagrams lose the meaning of quantum field theory objects and only display how quark 
flavours can be transferred form an initial hadron state to a final hadron state in hadron- 
hadron low-energy transitions. The coupling constants of such transitions desribed in 
terms of divergent parts of quark diagrams and depending of the cut-off A x = 0.94 GeV 
and the constituent quark mass m = 0.33 GeV can be expressed in terms of effective 
phenomenological coupling constants of low-energy hadron interactions given by effective 
chiral Lagrangians [6-9,17,23,24]. Hence, such a description of strong low-energy inter- 
action of hadrons can be valued as good established as the effective chiral Lagrangian 
approach. As has been shown in Refs. [6-10,14-16] (CHPT) g allows to develop consistent 
expansions in powers of both current quark masses and external momenta of interacting 
hadrons, i.e. chiral perturbation theory at both tree-meson and one-loop meson level, in 
good agreement with experimental data. 

For the matrix element < (vr + (g + )7r + (p + )7r _ (p_)) NR |n(0) 7m (1 — 7 5 ) d(0)|0 > we obtain 
the following 

zh/9 m 2 

< (n + (q + )n + (P+)*-(P-)hR\m 7 M (1 - 7 5 ) <*(0)|0 >= vt^ (p + + ?+)". (3.8) 

S 

This gives one 

2\/2 m 2 

p, < (7r + (g + )7r + (p + )7r-(^)) NR |n(0) - 7 5 ) d(0)|0 >= -5- (M 2 D+ + g 2 ),(3.9) 

^ S 

where we have set q\ = p\ = p 2 _ = and q 2 = (p + + q + ) 2 . 

Thus, in our approach the r.h.s. of Eq.( |3.9j ) is not proportional to the sum of current 
quark masses 

Pai < (7r + (q + )n + (p + )n-(p^ R \u(0) 7 "(1 - 7 5 ) d(0)|0 >= 

= -(m 0u + m od ) < (7r + (g + )7r + (p + )7T-(p_)) NR |M(0)7 5 rf(0)|0 >, (3.10) 

that can be expected when assuming the validity of the application of the equations of 
motion for the free current quark fields, i.e. u(0)p(l — 7 5 )rf(0) = — (m 0u + m d) , u(0)7 5 (i(0). 
This entails a vanishing of the matrix element p p < (7r + (g + )7r + (p + )7r~(]9_))NR|w(0)7 M (l — 
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7 5 )g?(0)|0 > in the chiral limit. Let us show now that at p 2 = M 2 + ^ M 2 the matrix 

s 

element p^ < (7r + (g + )7r + (p + )7r~(p_))NR|w(0) 7 M (1 — 7 5 )<i(0)|0 > does not vanish in the 
chiral limit. 

Let us consider a more general matrix element < n + (q + )n + (p + )n~(pJ)\u(<S)'y ll (\ — 
7 5 ) d(0)|0 > and represent it in terms of invariant amplitudes 

-i < n + (q+)n + (p+)n-(p_)\u(0) 7^(1 - 7 5 ) d(0)|0 >= 

= iq (p+ + g+f + F 2 p» + iF 3 e^pvp+aq+p, (3.11) 

where iq (i — 1, 2, 3) are invariant amplitudes free of kinematical sigularities. Multiplying 
both sides of Eq. ( |3.11|) by a momentum p^ we arrive at the expression 

-ip„ < n + (q + )n + (p + )n-(p^)\u(0)Y(l-l 5 )d(0)\0 >= 

= \Fi (p 2 - pi + q 2 ) + F 2 p 2 . (3.12) 

Imposing then the Adler condition, i.e. demanding a vanishing of the amplitude at p_ — > 
0, we obtain a relation between invariant amplitudes F 2 = —Fx, which ensues the equation 

- i lim p, < 7i + (q + )n + (p + )n-(p_)\u(0) 7 "(1 - 7 5 ) d(0)|0 >= i(iq + F 2 )p 2 = 0, (3.13) 

where we have used the relation q 2 —>■ p 2 valid in the limit p_ — > due to a conservation 
of 4-momentum p = g + + p + + p_. Using the relation i*2 = — iq we get 

- zp„ < 7r+(g + )7r + (p + )7r-(p_)|^(0) - 7 5 ) d(0)|0 >= ^ iq (g 2 - p 2 - p 2 _). (3.14) 

The correctness of this expression can be verified by matching Eq.( |3.14j) with the ampli- 
tude of the elastic 7r + 7r + -scattering (n + + tt + — ► n + + 7r + ). Setting p 2 = p 2 _ = M 2 and 
denoting g 2 = (p + + q + ) 2 = s we reduce the r.h.s. of Eq.( 3.14 ) to the form 

- ip, < n + (q + )n + (p + )n-(p^)\u(0) -f(l - 7 5 ) d(0)|0 >= ~ iq (s - 2 M 2 ), (3.15) 

which up to the common factor agrees well with the amplitude A4(ti + + 7r + — > 7r + + 7r + ) = 
(2 M 2 — s)/F 2 given by Weinberg. 

In the chiral limit p\ = q\= p 2 _ = M 2 —>■ we find 

< 7r + (g + )7r+(p + )7r-(p_)|n(0)7 /1 (l-7 5 )rf(0)|0 >= \f x s. (3.16) 

It is correct, since nothing suppresses the elastic 7r + 7r + -scattering for massless 7r + -mesons. 
This confirms the validity of Eq. (|3.9| ) and non-vanishing behaviour of the matrix element 
Pa < (7T + (g + )7r + (j9 + )7r-(p_)) NR |u(0)7^(l -7 5 ) d(0)\0 > in the chiral limit at p 2 = M 2 + ^ 
Ml 

Substituting Eq. (|3~9| ) in Eq. (|3.3|) we get the amplitude of the — > (ty + ty + 7y ) nr 
decay amplitude 



F n + Mi + + q 2 

M(Dj -> (tt+tt+tOnr) = d 2V2m 2 G F V* V ud — ^ . (3.17) 



The partial width of the D+ — > (ty + ty + 7y ) nr decay computed at M x = is given by 



T(Dj^ (n+n+n-) NR ) = = |Ci| 2 \G F V* s V ud \ 2 ~0- M D + 

= \d\ 2 x 0.87 x 10 10 s _1 , (3.18) 



where the numerical value has been obtained at m — 0.33 GeV, \V CS \ = 1.01, \V u d\ = 0.975 
and M n+ = 1.97 GeV. 

s 

Matching this value with the value of the partial width of the — > </>7r + decay we 
arrive at the ratio 

which agrees well with the experimental data Eq. ( |1 . 1| ) . 



4 Conclusion 

The application of the effective quark model with chiral U(3) x £7(3) symmetry incorporat- 
ing HQET and (CHPT)q to the computation of the partial widths of the — > 07r + and 
— > (7r + 7r + 7r~)NR decays has shown the agreement of the theoretical values with the 
experimental data. We have found T(D+ -> (7r + 7r+7r-) NR )/r(L)+ -> 0vr+) = 0.24 ± 0.12, 
which agrees well with the experimental value T(Z)+ — > (7r + 7r + 7r^)NR)/r(D+ — > 07r + ) = 
0.29 ±0.09 ±0.06. 

Thus, one can conclude that such an effective quark model describes reasonably well a 
low-energy dynamics of heavy-light meson interactions, and one does not need to include 
unnaturally heavy light current quarks in order to explain the experimental data on the 
7J+ -> (7T + 7r + 7r") NR decay [1]. 

Using the experimental probability of the Df — > 07r + decay, i.e. Br.(Df — > 07r + ) = 
(3.6 ± 0.9) %, we have estimated the time life of the D*-meson: t d + = (0.59 ± 0.15) x 
10 -12 s. This theoretical value agrees well with the experimental one: (r D +) exp = (0.467 ± 
0.017) x KT 12 s [11]. 

As has been discussed in Ref.[14] the theoretical uncertainty of our effective quark 
model is about 50%. However, for all cases of the application of this model, the resultant 
agreement between theoretical and experimental data turns out as usually much better. 

Recall that our results have been obtained for the factorized amplitudes. The step 
beyond the factorization approximation is to take into account one-meson loop contribu- 
tions. In (CHPT)g the procedure of the computation of one-meson loop corrections has 
been considered in Ref.[6] (Ivanov). We are planning to investigate such contributions in 
our forthcoming publications. 

5 Appendix. The 7r + — meson pole contribution 

In (CHPT)g the contribution of the 7r + -meson pole to the matrix element Eq.( |3.4| ) is 
given by [8] 

z < (n + (q + )n + (p + )n-(p_)) NR \u(0)^(l-^)d(0)\0> w+ ^ pole = 



Ng« fd 4 k J 5 1 5 1 1 1 ( Ngt\ fd A k 
: tr<7^7 - 7 - 



167T 2 J ix 2 i \ m -k+p m-k) P 2 - Ml\ IQtt 2 J J txH 

tr{ 7 5 :7 s J— - 7 5 1 -7 5 ^ T l + ( g+ ^ p+ ). (5.1) 

[ m — k + p m — k + p — p + m — k + q + m — k ) 

Taking into account that M D + > A x > we reduce the r.h.s. of Eq. ( |5.1| ) to the form 

lml 



i < (7r + (g + )7r+(p + )7r-(p_)) NR |M(0)7 /i (l-7 5 )^(0)|0> w+ _ pole = 



r d 4 k f „, 1 | 1 r d 4 k 
J ^ tr j 7 P ^kjp 2 -M 2 J ^ 



2V2 F n M 4 



D1 



ti'^7 5 1 7 5 1 



m — k + q+ m — k J 
+(q+~p+), (5.2) 

where we have used the relations g n = \/2m/F w , Ng 2 /%ii 2 = 1 and p 2 = M 2 + . 
The momentum integrals entering the r.h.s. of Eq. (|5.2p equal [6-9] 



fdk f 1 1 „ o „ 

/ — r-tr<7 M p =4ra »/, 

J nh [ m-k) 



d 4 k ( 5 1 5 1 1 

■tr<^7 5 —Y ->=4mv, (5.3) 



tt^ [ m — k + q + m — k 

where v = — < 0|gg|0 > / F 2 = 1.92 GeV [6-9]. This gives the following contribution of 
the 7r + -pole 



1 < (7r + (g + )7r + (p + )7r-(p_)) NR |M(0)7 /i (l-7 5 )rf(0)|0>^_ 



pole 



^i?4 p " (5 - 4) 

The contribution of the 7r + -meson pole to the amplitude Eq.( |3.3| ) is given by 

M(Dj - (7T+7r + 7r-) NR ) 7r+ _ polc = d 2^ m 2 G F l£ Kd -r^-. (5.5) 

Thus, we have shown that the contribution of the 7r + -meson pole is of order 0(1/M 4 +) 
compared with the contribution of the momentum integral Eq. (|3.5| ) and can be ignored, 
correspondingly. 
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